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Abstract 

Within finite-density QCD sum-rule approach we investigate the self-energies 
of E hyperons propagating in nuclear matter from a correlator of E inter- 
polating fields evaluated in the nuclear matter ground state. We find that 
the Lorentz vector self-energy of the E is similar to the nucleon vector self- 
energy. The magnitude of Lorentz scalar self-energy of the E is also close to 
the corresponding value for nucleon; however, this prediction is sensitive to 
the strangeness content of the nucleon and to the assumed density dependence 
of certain four-quark condensate. The scalar and vector self-energies tend to 
cancel, but not completely. The implications for the couplings of E to the 
scalar and vector mesons in nuclear matter and for the E spin-orbit force in 
a finite nucleus are discussed. 

PACS numbers: 24.85. +p, 21.65. +f, 21.80.+a, 11.55.Hx 

Typeset using REVTpX. 



*Permanent address: Instituto de Fi'sica, Universidade de Sao Paulo, 01498 - SP- Brazil. 



1 



I. INTRODUCTION 



Large and canceling isoscalar Lorentz scalar and vector self-energies for propagating nu- 
cleons in nuclear matter (unlike the conventional nonrelativistic picture) are essential for the 
successes of relativistic nuclear phenomenology ]!]-§]]. It has been shown recently, within a 
finite-density QCD sum-rule approach, that this physics might be motivated from quantum 
chromodynamics (QCD) |5|-^|. (Other applications of sum rule methods to finite density 
problems are discussed in Refs. JT0-O|). The predictions of QCD sum rules for the nucleon 
self-energies are found to be consistent with those obtained from relativistic phenomenologi- 
cal models (e.g., the relativistic optical potentials of Dirac phenomenology [|T],3 or Brueckner 
calculations [||,|J). However, it was found in Refs. J7||5| that the sum- rule predictions for the 
scalar self-energy are sensitive to assumptions made about the density dependence of certain 
four-quark condensates. Clearly, definite conclusions are not yet justified and further tests 
of the approach are needed. To this end, one is naturally led to consider the hyperons in 
nuclear matter within the same framework. 

By studying the hyperons, one can use both experimental data on hypernuclei and rela- 
tivistic phenomenology to confront with the QCD sum-rule predictions. Various investiga- 
tors |T6 H28|1 have applied the relativistic phenomenology to hypernuclear physics. In these 
relativistic models, the hyperons are coupled to the same scalar and vector fields as the 
nucleon, but with different coupling strengths. (These couplings are also of great relevance 
to other branch of physics |21| .) However, these strengths are not well-established. QCD 
sum-rule predictions for the scalar and vector self-energies may provide valuable insight into 
these couplings in the nuclear medium. 

In Ref. the self-energies of a A hyperon propagating in nuclear matter have been 
studied using the finite-density QCD sum-rule methods. The sum-rule calculations indicate 
that the self-energies of the A are only about ~ of the corresponding nucleon self-energies, 
suggesting a significant deviation from SU(3). In this paper, we further extend the finite- 



2 



density QCD sum-rule methods to explore the self-energies of X hyperons in an infinite 
nuclear matter. 

One of the compelling successes of relativistic models in describing nucleon-nucleus in- 
teractions is the naturally large spin-orbit force for nucleons in a finite nucleus. This force 
depends on the derivatives of the scalar and vector optical potentials, which add construc- 
tively. An analogous prediction for the A hyperon would seem to be problematic, if one 
adopts the naive SU(3) prediction that each coupling for the A should be | the coupling 
for the nucleon pl| , |32| , |20| -F25|1 . [That is, if one assumes that the scalar (a) and vector (u) 
mesons couple exclusively to the up and down quarks and not to the strange quark.] In 
the A-nucleus system, recent experiments indicate that the spin-orbit force is small, and 
even consistent with zero . This has raised questions about the validity of relativistic 

nuclear phenomenology for the hyperons. 

In Ref. [[TJJ, a weak A spin-orbit force was achieved by taking the potentials (i.e., the 
couplings) for the A to be a factor of three smaller than for the nucleon. More recently it 



has been suggested |17|j20| --p5| that larger values of the scalar and vector coupling strengths, 
consistent with the naive SU(3) prediction, can be used if a new tensor coupling between 
hyperons and the vector meson is introduced. In Refs. [p0| , p3| , ^4l| , it was argued that 
a quark-model picture implies that the tensor couplings of the hyperons (A,E,H) to the 
vector meson differ in their magnitudes and signs and hence their contribution to the spin- 
orbit force is different. For A, this picture leads to a tensor coupling with strength equal 
in magnitude to the corresponding vector coupling, but with the opposite sign. The net 
result, in combination with the scalar contribution, is a small spin-orbit force for the A. 
However, the experimental information available is not sufficient to resolve the effects of the 
tensor couplings. The sum-rule calculations in Ref. JHJ suggest that the coupling of the A 
to mesons are only about | of the corresponding nucleon couplings, implying a significant 
deviation from the naive SU(3) prediction and a weak A spin-orbit force; the tensor coupling, 
however, does not appear in the sum-rule calculations in a uniform nuclear matter. 
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At the moment, experimental evidence of E hypernuclei is insufficient. However, a num- 
ber of authors [E3-f26j have extended the relativistic phenomenology to E hypernuclei and 



presented theoretical predictions. In Refs. p3|-f25|, the naive SU(3) prediction for the cou- 
pling strengths of E to the scalar and vector mesons were adopted (the vector coupling for 
E is | the coupling for nucleon and the scalar coupling for E is slightly smaller than | the 
coupling for nucleon); in addition the tensor coupling between the E and the vector meson 



was included. In the quark-model picture of Refs. P^j3^j20|j24]1 , the tensor coupling of E to 
uj has the same sign as the corresponding vector coupling, in contrast to the A case. With 
the quark model values for the tensor coupling, the spin-orbit force for the E was found to be 
comparable with the nucleon spin-orbit force. In Ref. |29j the tensor coupling was omitted 
and universal couplings were assumed for all hyperons; the ratio of E to A spin-orbit force 
is about 0.9. The sum-rule approach may offer independent information on the scalar and 
vector couplings for E and new insight into deviations from SU(3) and into the spin-orbit 
force, although the prediction of a tensor coupling is not tested in the calculations described 
here. 

The finite-density QCD sum-rule approach focuses on a correlation function of interpo- 
lating fields, made up of quark fields, which carry the quantum numbers of the hadron of 
interest. The correlation function is evaluated in the ground state of nuclear matter instead 
of the QCD vacuum (as in the usual sum rules). For spin-| baryons, this function can, 
in general, be decomposed into three invariant functions of two kinematic invariants. The 
appearance of an additional invariant function compared to the vacuum case is due to an 
additional four-vector at finite density, the four-velocity of the nuclear medium. In the rest 
frame of the medium, the analytic properties of these invariant functions can be studied 
through a Lehmann representation in energy. 

The quasibaryon excitations (i.e., the quasiparticle excitations with baryon quantum 
numbers) are characterized by the discontinuities of the invariant functions across the real 
axis, which we use to identify the on-shell self-energies. By introducing a simple phenomeno- 
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logical ansatz for these spectral densities, we obtain a representation of the correlation 
function. On the other hand, the correlation function can be evaluated at large space- 
like momenta using an operator product expansion (OPE). By equating these two different 
representations using appropriately weighted integrals, one obtains QCD sum rules, which 
relate the baryon self-energies in the nuclear medium to QCD Lagrangian parameters and 
finite-density condensates f7|. 

We find that the finite-density QCD sum rules predict strong Lorentz scalar and vector 
self-energies for the E hyperon. In particular, the E vector self-energy is very similar to 
the sum-rule prediction for the nucleon vector self-energy. In terms of an effective theory 
of baryons and mesons, this implies larger couplings than would be indicated by SU(3) 
symmetry and the quark model. The predictions for the vector self-energy are essentially 
determined by the nucleon density and the fact (s^s) PJV = 0, and the predicted ratio of the 
vector self-energy to the zero-density E mass is found to be largely insensitive to the details 
of the calculation. In contrast, the predictions for the scalar self-energy are sensitive to the 
strangeness content of the nucleon and to the assumed density dependence of four-quark 
condensate (qq) 2 pN - If we assume that the nucleon has a significant strangeness content and 
the four-quark condensate (qq) 2 PN has weak density dependence, then the magnitude of the 
E scalar self-energy is found to be close to the corresponding value for nucleon and tends to 
cancel the vector self-energy, which is compatible with relativistic phenomenology. On the 
other hand, if the strangeness content of the nucleon is small or (qq) 2 PN depends significantly 
on the nucleon density, the magnitude of E scalar self-energy is very small and the net E 
self-energy is sizable and repulsive. 

This paper is organized as follows. The sum rules for E hyperons in nuclear matter are 
established in Sec. [0. The results are presented in Sec. [Ill] and discussions are given in 
Sec. [IV]. Section |VJ is a summary. 
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II. FINITE-DENSITY SUM RULES FOR S 



In this section, we derive the finite-density sum rules for E. We use the methods de- 
veloped in Refs. [|7|||, and refer the reader to these papers for further details. We work to 
leading order in perturbation theory. (The leading-logarithmic perturbative corrections are 
included through anomalous dimension factors.) In the OPE for E correlator, we consider 
all condensates up to dimension four, and the terms up to first order in the strange quark 
mass m s . Contributions proportional to the up and down current quark mass are neglected 
since they give numerically small contributions. In addition, we include the contributions 
from the dimension-six four-quark condensates (in the first order of m s , we only consider the 
terms which do not vanish in the zero-density limit). All other dimension six and higher- 
dimensional condensates are neglected since we expect their contributions to be small in the 
region of optimal Borel mass. 

QCD sum rules for E at finite density study the correlator defined by 

Ifc(ff) = i J d 4 x e ^(*o|T[^(^ E (0)]|*o> , (2.1) 

where f]T,(x) is a colorless interpolating field, constructed from quark fields, which carries 
the quantum numbers of E. The ground state of nuclear matter |^ ) is characterized by the 
nucleon density in the rest frame and the nuclear matter four- velocity u^; it is assumed 
to be invariant under parity and time reversal. Here we consider the baryon interpolating 
fields that contain no derivatives and couple to spin-| states only. There are two linearly 
independent fields with these features. For the nucleon, Ioffe [|35l , |36| has argued that the 
optimal choice is rj N (x)=e a t, c [u^(x)C'y tl Ub(x)] ■y 5 , y^d c (x) (for the proton). The interpolating 
field for E + can be directly obtained by an SU(3)-transformation of T)n(x) |J7J: 

Ve(x) = e abc [ul(x)C'j IM u b (x)} 7 5 7^s c (x), (2.2) 

where T denotes a transpose in Dirac space, C is the charge conjugation matrix, a, b, and 
c are color indices, and u(x) and s(x) are the up and strange quark fields, respectively. 
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The analogous interpolating field for S~ follows by changing the up quark field into down 
quark field. In the calculations to follow, we assume the isospin symmetry and use q to 
denote either up or down quark field. The interpolating field Eq. ( |2.2|) will be used in our 
calculations. 

Lorentz covariance, parity and time reversal then imply that the correlator can only have 
three distinct structures f||7]||: 

n E (g) = n s (g 2 , q-u) + U q {q 2 , q-u)i + U u {q 2 , q-u)i . (2.3) 

[A potential invariant function multiplying (q p u u — q u u p )a pu vanishes due to time- reversal 
invariance and a potential invariant function multiplying e xpfll/ (q\u p — q p u\)o pv vanishes due 
to parity conservation]. The three invariant functions, IT S , U q and U u , are functions of the two 
Lorentz scalars q 2 and qu. In the zero-density limit, U u — > 0, and II S and U q become functions 
of q 2 only. For convenience, we will work in the rest frame of nuclear matter hereafter, where 
u p = (1,0) and q-u — > qo; we also take Ui(q 2 ,q-u) — > TI^ (g , |q|) (i = {s, q, u}). To obtain 
QCD sum rules, we need to construct a phenomenological representation for Il s (g) and to 
evaluate Il s (g) using OPE techniques. 

In the rest frame of nuclear matter, the analytic properties of Hz(q) can be made manifest 



by a Lehmann representation [38], which leads to a dispersion relation in q of the form 



IT(go, |q|) = t^t I" ^ An ^' |q|) + polynomial (2.4) 
2m J-oo uj — q 



for each invariant function Hi, i = {s,q,u}. The polynomial stands for contributions from 
the contour at large |go|, which will be eliminated by a subsequent Borel transform (see 
below). The discontinuity AIT (which is the spectral density up to a constant) defined by 

AUi(u, |q|) ee lim ^(w + ie, |q|) - Ui(u - ie, |q|)] , (2.5) 

contains the spectral information on the quasiparticle, quasihole, and higher energy states. 

In QCD sum-rule applications, one parametrizes the spectral density with a small number 
of spectral parameters characterizing the resonances in the channel of interest (e.g., poles, 



residues, etc.). In vacuum, the spectral weights for baryon and antibaryon are related by 
charge conjugation symmetry and one usually parametrizes the spectral density as a single 
sharp pole representing the lowest resonance plus a smooth continuum representing higher 
mass states. At finite density, the ground state is no longer invariant under ordinary charge 
conjugation. Thus, the spectral densities for baryon and antibaryon are not simply related. 

The width of E in free space is small and can be ignored on hadronic scales. At finite 
density, the width of the E will be broadened due to strong conversions. Here we assume 
that the broadened width is relatively small on the hadronic scale and that a quasiparticle 
description of the E is reasonable. In the context of relativistic phenomenology, the E is 
assumed to couple to the same scalar and vector fields as the nucleons in the nuclear matter, 
and are treated as quasiparticles with real Lorentz scalar and vector self-energies. We follow 
Ref. and assume a pole approximation for the quasibaryon and take the spectral ansatz 
to be (higher-energy states are included in a continuum contribution, as discussed below) 

M* \* 2 r 

AU s (u 1 \q\) = -27ri^f 5(u-E q )-5(u-E q ) , (2.6) 

2E q l J 

\* 2 

AIL>, |q|) = -2™^ [S(u - E q ) - 6(u - E q )\ , (2.7) 
AH^, |q|) = +2«^f [6(u - E q ) - 8{u - E q )\ , (2.8) 



where A s is an overall residue. Here we have defined M£ = Mg + E s , E* = J M|. + q 2 , 



E q = T, v + y + q 2 , and E q = Y> v — y M| + q 2 , where Ms is the mass of S and S s and 
are identified as the scalar and vector self-energies of a E hyperon in nuclear matter. The 
positive- and negative-energy poles are at E q and E q , respectively. Since we want to focus 
on the positive-energy E pole, we approximate the spectral functions for the positive-energy 
E by a quasiparticle pole while suppressing contributions from the region of the negative- 
energy excitations. This is achieved by manipulating the parts of the correlator that are 
even and odd in q (see below). 

We now turn to the QCD expansion of n s (g), which is obtained by applying the oper- 
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ator product expansion to the time-ordered operator product in Eq. ( |2.1| ). In the present 
formalism, the correlator is studied in the limit that go becomes large and imaginary while 
|q| remains fixed (in the nuclear matter rest frame). This limit takes q 2 — > — oo with 
\q 2 /q-u\ — > oo, which satisfies the conditions discussed in Ref. |$9|] for a short distance ex- 
pansion. At finite density, the OPE for the invariant functions of spin-| baryon correlator 
takes the general form 0|| 

iW^-u) =Y.C l n (q\q-u)(d n ) PN , (2.9) 

n 

where (O n ) PN = (^olOnl^o)- The C % n (q 2 ,q-u) {i = {s,q,u}) are the Wilson coefficients, 
which depend on QCD Lagrangian parameters such as the quark masses and the strong 
coupling constant. The most important feature is that the local composite operators are 
defined such that all density dependence of the correlator resides in the condensates; the 
Wilson coefficients are then independent of density The operators O n are ordered 

by dimension (measured as a power of mass) and the C l n {q 2 ,q-u) for higher-dimensional 
operators fall off by corresponding powers of Q 2 = — q 2 . Therefore, for sufficiently large Q 2 , 
the operators of lowest dimension dominate, and the OPE can be truncated after a small 
number of lower- dimensional operators. 

The Wilson coefficients only depend on g M , and the ground-state expectation values of 
the operators are proportional to tensors constructed from the nuclear matter four-velocity 
the metric g^, and the antisymmetric tensor e KX ^ u . In Eq. ( [2.9|) we incorporate the 
contraction of g M (from the OPE) and (from the ground-state expectation values of the 
operators) into the definition of the Wilson coefficients C l n (q 2 ,q-u). Thus the dependence 
on q-u is solely in the form of polynomial factors. We have also suppressed the dependence 
on the normalization point [i. 

The Wilson coefficients can be calculated using the fixed-point gauge pE0| , pif and standard 



background-field techniques P^H^JHTJ. To obtain the Wilson coefficients one can apply 



Wick's theorem to the coordinate space time-ordered product in Eq. ( |2.1| ), retaining only 
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those contributions in which the quark fields are fully contracted, and using the quark 
propagators in the presence of the nonperturbative nuclear medium, given in Refs. |8|,|30l, 
for the contractions. 

Although the four-quark condensates have dimension six, their contribution to the baryon 
correlator are particularly important since the corresponding Wilson coefficients do not carry 
the large numerical suppression factors typically associated with loops j|7|]. By using the 
quark propagators given in Refs. |p|,|30| in the calculations, one includes the contributions 
from the four-quark condensates automatically with the condensates in their in-medium 
factorized form (See Sec. [ITI| for discussion of the factorization approximation). 

For convenience we split the invariant function Tli(qo, |q|) into two pieces that are even 
and odd in q : 

n 4 (g , |q|) = II? (q 2 , |q|) + g n°(g 2 , |q|) , (2.10) 
for i = {s, q, u}. The results of our calculations are 

n° = -^ 2 ) 2 H-i 2 ) + ^ 2 H-i 2 )(ss) PN - ^f 2 (^D q) PN 

n° = ^H~Q 2 ) {Wv) PN - (Js) PN ) - ^(ss) PN tfq) PN , (2.12) 
' (M-q 2 ) - A) (-\(u.Gf + (u-GY' 



144vr 2 V v ^ ' q 2 J \ tt 



PN 

2 



, 51n( q) 2—1 (ss) PN + — ^ </ , — 2 j w <-nv/>.x 



2 

PN 



-^{q ] q) PN ^s) PN , (2.13) 

n° = -^M-q 2 ) {(q ] q) P N + <*M,*) , (2.14) 



1 — 2 q 2 H~q 2 )(7(q^q) PN + (sh) PN ) , (2.15) 



12tt 
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n° = ^ln(-g 2 ) (m s (ss) PN - 16(qHD q) PN -A(sHD s) PN ) 

-^(<ffM; w + tfgW<aM w ) . (2-i6) 

where we have introduced the notation G 2 = G A U G A ^ U and G AftA = ^ KXfMU G A u , and have 
omitted the polynomials in g 2 and g 2 , which vanish under the Borel transform. Here we 
take {qOq) PN = ((uOu) PN + (dOd) PN )/2, with O a combination of Dirac matrices, gluon 
field tensors, and covariant derivatives. In Eqs. ( BTTTJ — ETTBH , we have neglected the terms 
proportional to (qiD q) PN since (qiD q) PN = \(m u + m d )pN — (see Ref. gj). 

In the present calculations, we have neglected all the terms proportional to the dimension- 
five quark and quark-gluon condensates. There is no difficulty in calculating the Wilson 
coefficients for these condensates. However, we have no reliable approach to determine 
these condensates. With the estimated values of the quark and quark-gluon condensates 
(containing only light quarks) given in Ref. |§, the dimension-five quark and quark-gluon 
condensates make numerically small contributions to the nucleon sum rules [f|. We expect 
that the contributions from the dimension-five quark and quark-gluon condensates are also 
small in the £ sum rules, so we neglect these terms. Contributions from three-gluon con- 
densates, which are dimension six, are also expected to be numerically small compared to 
those of four-quark condensates since the four-quark condensates enter at tree level whereas 
the three-gluon condensates involve loops 0. 

The QCD sum rules follow by equating the spectral representation of the correlator to 
the corresponding OPE representation. To improve the overlap of the two descriptions for 
QCD sum rules for hadrons in vacuum, one typically applies a Borel transform to both sides 
of the sum rules. For practical purposes in the present approach, the Borel transform can 
be applied using the operator B defined by 
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m<?, ioi)i - (-^) /w 2 . iii) - /(^ 2 . iii) . 

Q 2 = - 9o 2 , M 2 = ^. (2.17) 

n 

The result f(M 2 , |q|) depends on the Borel mass M. 
The finite-density QCD sum rules are then given by 

B[Uf (ql |q|) - E 9 n?(gg, |q|)]QCD = Bp? (ql |q|) - ^IT?^ | q |)] phcn , ( 2 .18) 

for i = {s,q,u}, where the left-hand side is obtained from the OPE and right-hand side 
from the phenomenological dispersion relations. To see how this manipulation suppresses 
contributions from the region of negative-energy excitations, we rewrite Eq. (|2.4j) as 

n*(gb, q ) = \ duu} — 2 — — + ^r~ du} — 2 — ~ > 2 - 19 

where we have omitted the polynomial term. We then obtain 

B[II? (g 2 , |q|) - (ql |q|)] P he„ = ± [°° du (u - E^AU^U, |q|) . (2.20) 

The negative-energy pole contribution is now suppressed by the factor (uj — E q ), which equals 
zero at u = E q . For Borel M near the S energy, contributions of higher-energy states to 
the integral are exponentially suppressed. 

Perturbative corrections ~ a™ can be taken into account in the leading logarithmic 
approximation through anomalous-dimension factor, L = H35| , where /i is the nor- 

malization point of the operator product expansion and Aqcd is the QCD scale parameter. 



In our numerical calculations, we take fi = 0.5 GeV and Aqcd — 0.1 GeV [|B . 

Finally, the contributions from higher-energy states are roughly approximated using the 
leading terms in the OPE, starting at an effective threshold; these contributions can be 
included by modifying the terms with positive powers of M 2 on the OPE side of each sum- 



rule equation as follows [35 
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M 2 -> M 2 E = M 2 (l — e~ s SM f 
M 4 -> M 2 E l = M 4 

M 6 -> M 6 E 2 = M 6 



1 _ e So/ M ' I _H L + 

\2M i M 2 



(2.21) 
(2.22) 

(2.23) 



where we define the continuum threshold Sq = u 2 — q 2 , with luq the energy at the contin- 
uum threshold. In principle, the effective thresholds are different for positive and negative 
energies and for the different sum rules. The former differences are critical in some sum rule 



formulations |p0|| , but are not numerically important in the present formulation. Further- 
more, the thresholds are relatively poorly determined by the sum rules and effects due to 
different thresholds in different sum rules may be absorbed by slight changes in the other 
parameters. In the present paper, we use a universal effective threshold for simplicity. 

With the spectral ansatz Eqs. (PUfl-flEgp and the OPE results of Eqs. (E1ID-( P11D , we 
obtain the finite-density sum rules for the S: 



\tM*e-W-^/ M2 



'"' M 6 E,L- 8/9 -^E 1 (ss) 



16vr 4 
m. — rr 



4vr 2 



Pn 



+ ^E q M 2 E ((Jq) PN - ( S t s)pjv L s/9 



2tt 2 

+ ^q 2 (^og) P ^- 8/9 + ^ 

- ^ 9 (^)p w (g t g) PJV , 
M 6 A/a M 2 



A * 2 e -(^-q 2 )/M 2 



2 

PAT 



(2.24) 



32tt 4 



+ 
+ 



M 2 



1447T 2 

m 



En 



327T 2 

M 2 



G 2 ) Eo^" 4/9 



71 



187T 2 

4M 2 
M 2 



+ 



9tt 2 
E„ 



PN 

(u-G) 2 + (u-Gf 

(ss)p N L^ 9 

(shD s) PN L-^ 



£-4/9 



PN 



M 5E> 



' M 2 



E 



M 2 



M 2 E ((Jq) PN + (sh) PN ) 



6n 2 

+ lmi N L 4/9 + - 3 ( q k)p N (s ] s) PN L-^. 



[2.25) 
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- §^M 2 E (m s (ss) PN - 16(qUD q) PN - 4(sUD s) PN ) L 

/^[( m . G )2 + (m .g) 2 1\ l- 4 / 9 



36tt 

+ ^Vff>, w ((Q ] l) PN + (s ] s) PN ) . (2.26) 

Here E , Ei and E 2 have been denned in Eqs. ( |2.21| )-( pi.23| ). We have ignored the anomalous 
dimensions of dimension four operators, either because the operators are renormalization- 
group invariant, because the anomalous dimension is small, because the corresponding con- 
densates give small contributions, or because the accuracy to which the nucleon matrix 
elements of the operators are known is such that anomalous- dimension corrections represent 
an unwarrantable refinement (see Sec. pT|) . We will adopt the values of the corresponding 
condensates at the scale of 1 GeV in our numerical calculations. The four-quark opera- 
tors are in general not renormalization-covariant, so they mix with one another under the 
renormalization group [37]. In vacuum, the anomalous dimension effects do not violate the 



factorization assumption to within 10% [J47| , and thus one usually assumes that the anoma- 
lous dimension of a four-quark operator is equal to the sum of the anomalous dimensions of 
the factorized operators |^5|,^TJ. In the present paper, we follow this assumption. 

III. RESULTS 

A. In-medium condensates 

To obtain the predictions for the E self-energies from the finite-density sum rules derived 
in the previous section, we need to know the in-medium condensates appearing in the sum 
rules. To first order in the nucleon density, one can write 

(d) PN = (6) vac + (d) NPN + ■ • ■ , (3.1) 
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where • • • denotes correction terms that are of higher order in pjy, and (0)n is the spin 
averaged nucleon matrix element. (Note that this is not a Taylor series in pjy.) For a 
general operator O there is not a systematic way to study contributions to (0) PN that 
are of higher order in pjv- Model-dependent estimates in Ref. [EJ suggest that the linear 
approximation to (qq) PN should be good (higher-order corrections ~ 20% of the linear term) 
up to nuclear matter saturation density. Here we assume the first-order approximation of 
all condensates to be reasonable for calculating scalar and vector self-energies up to nuclear 
matter saturation density. Justifying the limits of this type of density expansion is an 
important topic for further study. 

The simplest in-medium condensates are {q [ q) PN and (s^s) PN . Since the baryon current is 
conserved, (q^q) PN is proportional to the (net) nucleon and strangeness densities: (q^q) PN = 
|pjv, and (s' ! s) PN = 0. These are exact results. The other dimension three and four quark 
and gluon condensates have been studied in Refs. ||10|||30|]. The results are: 



(qq) Pn = (qq)^ + (qq)NpN = (gg) vac + ^pn , (3.2) 

(ss) PN = (ss) vac + (ss) N p N = (ss) v3jC + y-^^PN , (3.3) 
^1 G A =/^l G 2 ) -(650MeV) PAr , (3.4) 

7T / pN \ 7T /vac 

— \(u'-G) 2 + (u'-G) 2 }) = -(lOOMeV)pjv , (3.5) 

7T L J / PN 

(qUD q) PN = (180MeV) P7V , (3.6) 
(sUD s) PN = f(ss) PN + (18MeV) PiV . (3.7) 

where m q = \{m u + ma) is the average of the up and down current quark masses, cr^v is 
the nucleon a term, and y = (ss) n/ (qq) n is a real parameter, measuring the strangeness 
content of the nucleon. 

Four-quark condensates are numerically important in both the vacuum and the finite- 
density baryon sum rules because they contribute in tree diagrams and do not carry the 
numerical suppression factors typically associated with loops. In the sum rules derived in 
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Sec. U, we included the contributions from the four-quark condensates in their in-medium 
factorized forms; however, the factorization approximation may not be justified in nuclear 
matter [[^-|9|j30|| . In the nucleon sum rules, the "scalar-scalar" four-quark condensate (qq) 2 



PN 



gives important contributions; the sum-rule predictions are sensitive to the value of this 
condensate 0,3- This four-quark condensate also appears in the £ sum rules. In its fac- 
torized form, this scalar-scalar four-quark condensate has a strong density dependence; one 
might suspect that this strong density dependence is an artifact of the factorization approx- 
imation. Thus we follow Ref. |9j and parametrize the scalar-scalar four-quark condensate 
so that it interpolates between its factorized form in free space and its factorized form in 
nuclear matter: 

m 2 PN — m N = (i - fm)L + mr PN , (3.8) 

where / is a real parameter. The density dependence of the scalar-scalar four-quark con- 
densate is thus parametrized by /, and the density dependence of (qq) PN [see Eqs. (|3.2|) 1. 
The factorized four-quark condensate (qq)p N appearing in Eqs. (|2.24| - pl26D will be replaced 
by (qq) 2 PN in the calculations to follow. Studies of nucleon self-energies with QCD sum rules 
yield results in strong contradiction to experiment unless < / < 0.5 [0,0; however, there is 
not an independent determination of the density dependence. Here we will consider only this 
range of /. The other four-quark condensates have much smaller numerical contributions, 
so we just use the in- medium factorized form for simplicity. 

The values of and y remain controversial [p2|-[5^]. In this paper, we take a at = 45 MeV, 
which is the value obtained in a recent analysis ||54|| , and consider values of y in the range of 



— 0.6, which covers the values discussed in |p^-|5^]. For the condensates in vacuum, we use 
(gg) vac ~ -(245 MeV) 3 (m q ~ 5.5 MeV) |7]||, and take (ss) vac = 0.8(gg) vac and 



((a s /7r)G 2 ) vac = (330 MeV) 4 [P7yT4"f . We consider the strange quark mass m s in the range 
of 100 - 200 MeV. 
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B. Sum-rule analysis 



In principle, the predictions based on the sum rules should be independent of the aux- 
iliary parameter M 2 . In practice, however, we have to truncate the OPE and use a simple 
phenomenological ansatz for the spectral density; thus one expects the two descriptions to 
overlap only in some limited range of M 2 (at best). As a result, one expects to see a "plateau" 
in the predicted quantities as functions of M 2 . The studies of the vacuum sum rules for the 
octet baryons and the nucleon sum rules at finite-density show that the sum rules truncated 
at dimension-six condensates do not provide a particularly convincing plateau H^fl.Q 
Nevertheless, we will assume that the sum rule actually has a region of overlap, although 
imperfect. We follow Refs. ||7|fl|30| and rely on the cancellation of systematic discrepancies 
by normalizing all finite-density self-energies to the zero-density prediction for the mass. 
One hopes that this might compensate for general limitations of the sum rules. All the 
finite-density results presented are obtained at nuclear matter saturation density, which is 
taken to be p N = (110 Me V) 3 . 

To analyze the sum rules and extract the self-energies, we sample the sum rules in the 
fiducial region of M 2 , where the contributions from the highest-order condensates included 
in the sum rule are small and the continuum contribution is controllable. Here we choose 
1.0 < M 2 < 1.6 GeV 2 as the optimization region [the ratios of the self-energies to the S 
mass are insensitive to the choice of the upper bound of the optimal region (see Fig. |6|)]. 
The study of the £ sum rules in vacuum suggests that the sum rules are valid in this 
region |)T],[55|. To quantify the fit of the left- and right-hand sides, we use the logarithmic 



measure 



"maximum^e-^ 2 -^ 2 ^ 2 , IE S /M*, IT, U'JE V } 



5(M 2 ) = In 



minimum{A^ 2 e-( £2 - q2 )/ M2 , WjM^, IE, U'JE V } 



(3.9) 



including direct-instanton effects in nucleon sum rules in vacuum leads to a more convincing 
plateau 

17 



which is averaged over 150 points evenly spaced within the fiducial region of M 2 . Here U' s , U' q , 
and denote the right-hand sides of Eqs. Q2.24|) -( p.26j ), respectively. The predictions for 



Mg, E„, uj 2 , and A s 2 are obtained by minimizing the averaged measure 5. To get a prediction 
for the E mass in vacuum, we apply the same procedure to the sum rules evaluated in the 
zero-density limit. 

In Fig. m, we displayed the optimized results for the ratios M^/M^ and E^/M^ as func- 
tions of y for m s = 150 MeV, |q| = 270 MeV, and three different values of /. (The momen- 
tum dependence of the self-energies for momenta below the Fermi surface is very weak, see 
Fig. |^). It can be seen that E„/M s is insensitive to both y and /. However, the M^/Ms 
varies rapidly with y and /. Therefore, the sum rule prediction for the scalar self-energy is 
strongly dependent on the strangeness content of the nucleon and on the density dependence 
of the four-quark condensate. For / = and the values of y in the range 0.4 < y < 0.6, the 
predictions are: 

M s /M s ~ 0.78-0.85 , (3.10) 
E„/M E ~ 0.18-0.19. (3.11) 

On the other hand, for / = and small values of y (0 < y < 0.2), we find E^/Ms ~ 0.18 
and M^/Ms ~ 0.92-0.98. As / increases, M^/Ms increases, which implies an even smaller 
magnitude of the scalar self-energy. The predictions for the ratios A^ 2 /A| and Sq/s also 
depend on y and /. For / = and large values of y (0.4 < y < 0.6), both the continuum 

2 

threshold and the residue A s are close to their corresponding vacuum values. For / = and 
small values of y (0 < y < 0.2), the continuum threshold increases by about 20% relative to 
the vacuum value and the residue A s 2 increases by about 50% relative to its vacuum value. 
As / increases, both the continuum threshold and the residue increase. 

One can see, from the sum rules Eqs. ( |2.24| )-( p.26[ ), that the ratios H^/Il' and n^/n' ? 



give and E„ as functions of Borel M 2 , and n' s /n^ in the zero-density limit yields Ms as 
a function of M 2 . In Fig. 0, the ratios M^/Ms and E,;/M s are plotted as functions of M 2 



for y = 0.5 and various values of /, with E q , E q and the continuum threshold fixed at their 
optimized values. The curves for the ratios M^/Ms and £„/M£ are flat, and thus imply a 
weak dependence of the predicted ratios on M 2 (though the individual sum-rule predictions 
before taking ratios are not flat). 

We plot A^ 2 e- ( ^ 2 - q2)/M2 , IT S (M 2 )/M£, Tl' q (M 2 ), and U' U (M 2 )/E V as functions of M 2 for 
y = 0.5 and / = in Fig. 0(a), with the predicted values for M£, £„, Uq, and A^. If the 
sum rules work well, one should expect to see that the four curves coincide with each other. 
It is seen that their M 2 dependence in the Borel region of interest turns out to be equal up 
to 15%. The overlap of the corresponding vacuum sum rules (i.e., the zero-density limit) is 
illustrated in Fig. 0(b). We observe that the quality of the overlap for the finite-density sum 
rules is similar to that of the corresponding sum rules in vacuum. 

The sensitivity of our predictions to m s is illustrated in Fig. [|, where y and |q| are 
fixed at 0.5 and 270 MeV, respectively. The predictions for M^/M-^ and T, v /M-^ are largely 
insensitive to changes in m a , with a variation of S^/M^ less than 10% in the range of 
m s = 0.1—0.2 GeV. Fig. [| shows the three-momentum |q| dependence of the predicted ratios. 
The results are only weakly dependent on three momentum in the range of |q| = 0—500 MeV. 
Finally, the dependence of the results on the choice of the upper bound of the Borel window 
is shown in Fig. [!]. The two dashed curves are obtained using a fixed Borel window at 
1.0 GeV 2 < M 2 < 1.6 GeV 2 . The two solid curves are the results obtained by taking an upper 
bound of the Borel window such that the continuum contributions to the phenomenological 
sides do not exceed 50% of the total phenomenological contributions to the sum rules (i.e., 
the sum of the quasiparticle pole and the continuum contributions) while fixing the lower 
bound at 1.0 GeV 2 . We have used the same procedure in extracting the £ mass in vacuum. 
We see that the two set of curves are almost indistinguishable, implying that changing the 
upper limit of the optimum Borel region does not affect the sum-rule predictions for the two 
ratios. 
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IV. DISCUSSION 



We note that the sum-rule predictions for the scalar self-energy are quite sensitive to 
the strangeness content of the nucleon and to the undetermined density dependence of 
certain four-quark condensate. Therefore, we cannot draw definite conclusions about the 
E scalar self-energy at this point. Nevertheless, we emphasize that the sum-rule prediction 
for the normalized vector self-energy, E„/Mx, is apparently insensitive to the details of 
calculations. For typical values of the relevant condensates and other input parameters, 
E v /M s ~ 0.18-0.21. The finite-density nucleon sum rules predict K/M N ~ 0.25-0.30 
Thus, we find (E„) s /(E 1 ,) A r ~ 0.8 — 1.1. 

This result, if interpreted in terms of a relativistic hadronic model, would imply that the 
coupling of the E to the Lorentz vector field is very similar to the corresponding nucleon 
coupling in the same ratio. This compares to the naive SU(3) prediction of |, which is 



obtained by assuming that the mesons couple directly to constituent quarks |31] , |20| -p2]] , and 
thus suggests a significant deviation from SU(3) symmetry in nuclear matter. This deviation 
can be attributed to two sources. First, the nuclear matter ground state is not SU(3) 
symmetric due to the absence of net strangeness. This leads to further deviation of various 
condensates from SU(3) relative to the deviation in vacuum. The second source originates 
from the baryon interpolating fields. In the present work, we used the interpolating field 
defined in Eq. ( |2.2| ), corresponding to an axial vector diquark, composed of two up quarks, 
coupled to a strange quark. However, these quarks contribute to the sum rules differently. 
This can be seen from the leading-order terms in the sum rules. In Eq. ( |2.26| ) there is an 
extra factor of seven multiplying (q^q) PN relative to the term proportional to (s^s) Pjv . Since 
the strange quark does not couple to the nuclear vector current (i.e., (s^s) PN = 0), the 
leading-order term of Eq. ( |2.26| ), which sets the scale for the E vector self-energy, is very 
close to the corresponding leading term for the nucleon. 

The dependence of the E scalar self-energy on y comes mainly from the leading-order 
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term (proportional to (ss) PN ) of Eq. Q2.24|) and the parametrization Eq. ( |3.3| ). If we assume 
that the nucleon has a large strangeness content (i.e., y ~ 0.4 — 0.6) and the four-quark 
condensate (qq) 2 PN depends only weakly on the nucleon density (i.e, if / ~ 0), we find 
M^/M E ~ 0.77 - 0.84, which implies S S /M E ~ -(0.16 - 0.23). With the nucleon sum- 
rule prediction M* N /M N ~ 0.65 - 0.70 0J|, we obtain (£ a ) E /(£ a )jv ~ 0.6 - 1.0. In a 
hadronic model, this implies again a coupling of the S to the Lorentz scalar field close to 
that for the nucleon. In this case, there is a significant degree of cancellation between the 
scalar and vector self-energies, which is compatible with that implemented in the relativistic 
phenomenological models. 



In contrast, if the strangeness content of the nucleon is small (i.e., y < 0.2) or if (qq) 2 



PN 



has a significant dependence on the nucleon density, the predicted ratio M E /M E is close to 
unity, implying that the scalar self-energy is very small. The predicted vector self-energy, 
on the other hand, is still essentially the same as the nucleon vector self-energy. Thus, in 
this case the sum rules predict incomplete cancellation, and hence a sizable repulsive net 
self-energy for the E. This result contradicts with that of the relativistic models. 

We now turn to the £ spin-orbit force in a finite nucleus. In the present approach, all sum- 
rule predictions are obtained for uniform infinite nuclear matter; thus one cannot obtain a 
direct prediction for the £ spin-orbit force. However, we can still get some partial information 
on the £ spin-orbit force by adopting an approach of Dirac phenomenology . In this 
approach, the coordinate-space potentials entering the Dirac equation for a baryon scattering 
from a finite nucleus are assumed to follow a Fermi distribution with two overall potential 
depths (scalar and vector), which are independent of nuclei, and which can be associated 
with the self-energies in infinite nuclear matter. The spin-orbit force follows by recasting the 
Dirac equation in Schrodinger form. The resulting spin-orbit potential is proportional to the 
sum of the magnitudes of the potential depths (that is, the scalar and vector self-energies) 
multiplied by the derivative of the assumed Fermi distribution [[J. 

The sum-rule predictions for the £ scalar and vector self-energies in the present paper 
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imply that the E spin-orbit force in a nucleus is somewhat weaker than (but comparable 
with) that felt by a nucleon, but much stronger than that felt by a A. This is consistent with 



that obtained in Refs. fl23|-[25|j. However, we emphasize that in Refs. p3^p5[ the scalar and 
vector couplings, consistent with SU(3), have been adopted and it is the extra tensor coupling 
of £ to the vector meson that enhances the spin-orbit force. Our sum-rule calculations, on 
the other hand, suggest that it is the strong scalar and vector couplings, deviated from the 
naive SU(3) prediction, that lead to large spin-orbit force. We also note that our sum- rule 
predictions in this paper and those in Ref. JJIJ for A do not agree with the universal coupling 
assumption (i.e., all hyperons couple to the scalar and vector fields with the same strength) 
suggested in Ref. p6|. 



V. SUMMARY 

In this paper, we have applied finite-density QCD sum-rule methods to investigate the 
self-energies of a E hyperon in nuclear matter. The approach focuses on a correlator of £ 
interpolating fields, evaluated in the nuclear matter ground state. A QCD expansion of the 
correlator is obtained by applying the operator product expansion. We retained the con- 
tributions from all condensates up through dimension four and to first-order in the strange 
quark mass m s , and we also included the contributions from the four-quark condensates. 
The expansion requires new condensates not present at zero density, as well as informa- 
tion on the density dependence of condensates. In the rest frame of the nuclear matter, a 
Lehmann representation for the correlator, with fixed three-momentum, leads to a dispersion 
relation in the energy variable. A simple quasiparticle pole ansatz, with real Lorentz scalar 
and vector self-energies, is assumed for the spectral functions associated with the three in- 
variant functions comprising the £ correlator. Contributions from higher-energy states are 
roughly approximated by a perturbative evaluation of the correlator, starting at an effective 
threshold. 
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The sum-rule analysis indicates that the S vector self-energy is similar to the corre- 
sponding nucleon self-energy. If we interpret this result in terms of a relativistic hadronic 
model, it implies that the vector coupling for the S is similar to the corresponding nucleon 
coupling to the vector meson. The vector self-energy is largely insensitive to the details of 
calculations, and is essentially determined by the nuclear matter density and the fact that 
(^s) PN = 0. 

The sum-rule predictions for the S vector self-energy, along with that for the A in 
Ref. PC|| , indicate a significant deviation from SU(3) in nuclear matter. This deviation 



is mainly due to the violation of SU(3) in the baryon interpolating fields and due to the 
violation of SU(3) in the nuclear matter ground state. At this stage, it is unclear whether 
the violation of SU(3) in the baryon interpolating fields is connected to the violation of 
SU(3) in the baryon wave functions. Further study of the dependence of the deviation on 
the choice of the interpolating field will be important in understanding the deviation from 
SU(3). 

The scalar self-energy is found to be quite sensitive to the strangeness content of the 
nucleon and the assumed density dependence of the four-quark condensate (qq) 2 PN - We 
parametrized the density dependence of this four-quark condensate in terms of its factorized 
form in free space and its factorized form in nuclear matter. If the nucleon has a significant 
strangeness content and if the four-quark condensate (qq) 2 PN has a weak density dependence, 
the sum rules predict strong scalar and vector self-energies. This is qualitatively compatible 
with relativistic models. On the other hand, if the strangeness content of the nucleon is 
small or ^qq) 2 PN depends significantly on the nucleon density, the sum rules predict a very 
small scalar self-energy and a strong vector self-energy, which differs from the relativistic 
models. Clearly, further study of the strangeness content of the nucleon and the four-quark 
condensates in the nuclear medium is very important, along with analyses of the higher- 
order density dependence of other condensates and the contributions from the condensates 
with higher dimension. 
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The sum-rule predictions for the £ scalar and vector self-energies seem to imply a some- 
what weaker spin-orbit force for £ in a nucleus than that for a nucleon, but much stronger 
than that for a A. This is compatible with those obtained from relativistic phenomenologi- 
cal models with an extra tensor coupling between the hyperons and the vector meson. This 
possible tensor contribution, however, cannot be observed in uniform nuclear matter. 

It is a straightforward exercise to study the self-energies of S within the same framework. 
However, in the S sum rules, the leading order term of Il u is small relative to the higher 
order terms. Thus, it might be difficult to extract useful information since one does not 
have much control over the values of the higher order terms. 
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FIGURES 

FIG. 1. Optimized sum-rule predictions for M^/My, and Y, v /M-£ as functions of y. The three 
curves correspond to / = (solid), / = 0.25 (dashed), and / = 0.5 (dotted). The other input 
parameters are described in the text. 

FIG. 2. Ratios M|/Ms and £„/M s as functions of Borel M 2 for y = 0.5, with optimized 
predictions for E q , E q , and the continuum thresholds. The three curves correspond to / = 
(solid), / = 0.25 (long-dashed), and / = 0.5 (dotted). The other input parameters are the same 
as in Fig. [I]. 

FIG. 3. (a) The left- and right-hand sides of the sum rules as functions of Borel M 2 for y = 0.5 

2 

and / = 0, with the optimized values for M£, T, v , Sq, and A^ . The other parameters are the same 
as m Fig. |l|. The four curves correspond to H' s /M^ (solid), (dashed), U' u /T, v (dot-dashed), and 
A^e - ^ -0 ! 2 )/^ (dotted), (b) The left- and right-hand sides of the corresponding vacuum sum 
rules. The three curves correspond to U' s /M-£ (solid), 11^ (dashed), and A|e _M s/ M2 (dot-dashed) 
at the zero-density limit, with the optimized values for Ms, so, and Af;. 

FIG. 4. Optimized predictions for M^/Ms and £„/Ms as functions of m s for y = 0.5. The 
three curves correspond to / = (solid), / = 0.25 (dashed), and / = 0.5 (dotted). The other input 
parameters are the same as in Fig. |l[ 

FIG. 5. Three-momentum dependence of the predicted M^/Mx and T, v /Ms for y = 0.5. The 
three curves correspond to / = (solid), / = 0.25 (dashed), and / = 0.5 (dotted). The other input 
parameters are the same as in Fig. |l[ 

FIG. 6. Optimized sum-rule predictions for M|./Ms and £„/Ms as functions of /, with y = 0.5. 
The other input parameters are the same as in Fig. 0. The solid curves correspond to the results 
obtained by requiring the continuum contributions to be less than 50% in the fiducial Borel re- 
gion and the dashed curves correspond to the results obtained using a fixed Borel window at 
1.0 GeV 2 < M 2 < 1.6 GeV 2 . 
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